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1. Introduction 

It is well known that in lattice gauge theory the vertices are quite complicated especially for 
the pure gluon sector. This is because, on the lattice, we are preserving gauge invariance at finite 
cutoff and the lattice itself breaks the Lorentz symmetry explicitly. This is the main difficulty of 
lattice perturbation theory. To reduce the risk of errors and to alleviate the tedious task of deriving 
the vertices, it is desirable to have an automatic method. A first attempt was made by Luscher 
and Weisz about twenty yeas ago [^. They worked in momentum space and performed some 
calculations by using their algorithm which was restricted to closed loops sufficient for pure gauge 
theory. 

Recently a new algorithm, that we call bottom up algorithm, was proposed by Hart et. al.lQ]. 
A crucial point in this generalization is that it can deal with any parallel transporter, not only 
with closed loops. This allows to also include a fermion action or even a smeared HQET action. 
A relevant assumption for this algorithm has been translation invariance. Our main concern in 
this note is to extend the bottom up algorithm to the Schrodinger functional (SF)[Q] where this 
invariance is broken for the time direction. Before going to the extension we first summarize the 
position space version of the algorithm on the usual translation invariant lattice in the next section. 

2. Bottom up Algorithm 

In Ref.||^, the authors explain the algorithm in momentum space but here we will move to 
coordinate space. In the following we denote the antihermitean gauge fluctuation field by q^{x), 
and the link variable (still without background field) hy U {x,pi) = exp(^()^^(jc)). 

A first important point for the automatic operation is how to represent the vertices in a program. 
We consider the parallel transporter along a curve ^ on the lattice. The r'th order coefficient in 
the Taylor expansion Pr of the parallel transporter P[^;q\ is written as 

P[^;q] = l^^Pr[^;q], Pr= I (2.1) 

r=0 au...,ar k=\ 

where a = {pi,x) is a combined index labelling links, and a is color. We call "Toa^-ar = ^ai " ' ' ^o, ^ 
color factor, and fk an amplitude. The latter corresponds to a value of a reduced vertex with index 
{al , • • • , a[). The reduced vertex here means a vertex without the color factor. The information 
about the reduced vertex of order r is now encoded in a list consisting of Nr lines 

L^l^ = {al---,al;fk), {1.1) 

which holds the index structures and amplitudes of the corresponding reduced vertex. The list 
/.(O = {L^j^^k = \,..,Nr} corresponds to the all nonzero elements of the reduced vertex of order 
r. It will be constructed recursively by the multiplication algorithm to be given below. Note that 
in eq.( |2.1| ) the color factor is located outside of the ^-summation, because it is independent of the 
shape of the transporter. However, for the SF this nice structure is rendered more complicated by 
the background field as we will see in the next section. 
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As an explicit example for a list, let us look at the case of the single gauge link variable. From 
the expanded form we can obtain a hst L^''^ of the link variable, 

f/(^,M) = L^ L '^«i-a.<---<, lW = ( ;1), a = {^,x). (2.3) 

r elements 

It consists of only one line {Nr = 1) and is the elementary building block of the algorithm, that is, 
something like an initial condition. 

So far we have defined the fundamental elements on which the algorithm operates. Next we 
consider the case that a parallel transporter P is composed from P' and P" , P = P'P" and want to 
obtain the vertices of P from those of P' and P" . In other words, the problem is how to obtain a set 
of lists of P (up to a certain order), from those of P' and P" , {L/}p/ x {Lj}pii — > {Lk]p- Since any 
parallel transporter is composed of elementary one-link variables, by repeating the procedure, we 
can obtain the hsts for arbitrary parallel transporters. This is the origin of the name 'bottom up'. 
The algorithm can be easily understood by looking at an actual multipUcation of the coefficients of 
the Taylor expansion. The coefficient of P with order r, Pr, is expressed by those of P' and P" as 

p - y - p'p" 

r N' N"_ 

^a\...ar |/ \| 2^ ^ai -"^aj ^Oj+i -"^ar JiJj 
ai...ar s=0^-y 1=1 ;=1 

= ^ '^ai...ar ^ (la\ ' ' ' ^a*/^:- 
ai...ar k=l 

After the second equal-sign we inserted the exphcit form of the coefficients, and use the fact that 
the color factor is independent of the shape of the parallel transporter In the last step we combine 
the three summations into that over k. More precisely, we did a relabelling of the indices, and 
rewrote the amplitude factor. Finally the resulting list of P, is created by putting the new label 
structure and the new amplitude. The algorithm is summarized as 

• Relabehng: {«/,•••, «f , «] > ' • • , cc'f} — > {«/,•••, oct\ «0 

• Amphtude part: jf^pr^flfj — > fk 

• Creating hst: L[ x L'j — ^Lk = {al,---, a[;fk) 

This procedure should be carried out for < j < r, 1 < j < A'^' and I <j < N'/_^ if order r is desired. 
The algorithm has been implemented in the python script language, which is good at dealing with 
the complicated list operation. In this way, one can obtain the vertices for any parallel transporter. 

3. Extension to the Schrodinger Functional 

The essential new ingredient in the SF is the presence of the non- vanishing back ground field. 
It is induced by non trivial Dirichlet boundary conditions in the time direction. In this case, the hnk 
variable is expressed by the background field and the fluctuation field q^, 

U{x,^)=V{x,n)e^''"^'^''\ (3.1) 
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The main difficulty due to the presence of the back ground field occurs in the color factor. In the 
color factor, the back ground field can appear "randomly" between SU(3) generators, for example 

^a,a2a3-- = TaiTa2Taj- ■ ■ ^Ta^VTa^V ^Ta^---, On SF, (3.2) 

where arguments of Vs and powers {V or V^^) depend on the location of the links and its orien- 
tation. At first glance, it seems difficult to deal with it in a systematic way. Furthermore, since the 
background field depends on the Lorentz and position indices, the color factor is located inside of 
the sum over / and does not factorize anymore, 

Pr= I '^a,-a.l,qai---qT^fk^ £ I < ...«>f • • • ^.f A, on SF. (3.3) 

fll,...,a,. k=l a\,...,ark=\ 

One can solve the problem however for the restricted class of background fields that have 
mainly been used in applications of the SF [^], V{x,0) = 1, V{x,k) ='V{xq). In addition, we take 
the background field to be abelian, given by diagonal color matrices. This means that the generators 
in the V{xq) are written as linear combinations of the elements of the Cartan sub-algebra Hj, that 
is, V{xq) = e''^j^J^^°^^' with coefficients hj{xo). Now we find the nice equation, V{xo)IaV^^{xQ) = 
l^g^^a{xo) ^ and the la (Cartan basis) do not mix with other basis elements. This is because, in the 
adjoint representation, an element of the Cartan basis is an eigenstate of the Cartan generator, 
[Hj,Ia\ = pijala- The eigenvalue is a root, and the phase is a linear combination of the roots and 
hj{xo), (pai^o) = L; l^-jahjixo). Another property of the standard background fields that is relevant 
here is that the back ground field and the phase has a simple dependence on time, V{xQ + At) = 
y(.)co)exp (/Af(^), (pai^o +Af) = (pai^o) +A?i//o, where is the color electric field. This however 
needs to be relaxed later, see in next section. 

Our main finding is that by making use of the above properties of the background field any 
color factor of order r can be written as 

<....,. (^o) = [4, • • ■Iay^^{xo)e''^'^] ei^M^y^.M^l:\ (3.4) 

: ' {7(1) phase factor 
3x3 matrix 

where the background field has been moved to the right in the 3x3 matrix part. Actually we 
can show by induction that A and B are single component integer, and C and D are r component 
integer valued vectors. In the expression, the information of the lattice size and the back ground 
field are encoded in V, S, <p and Y- The lists are independent of these values, only xo,A,B,C,D 
are required. The former are only needed when implementing the vertex in a diagram calculation 
program at a second stage. Note that we need xq, since there is no translation symmetry for the 
time direction. We choose xq as the time component of the position of the left most link variable in 
the parallel transporter. The benefit of the expression is that we can separate the information about 
the list (A,B,C,D) and the lattice size and the back ground field. Therefore, we can do a symbolic 
list operation, independently of the details of the lattice and the background field. 

We have obtained a manageable expression for the color factor. Next we formulate the multi- 
plication for this structure of integer lists {xo,A,B,C,D). From an actual multiplication of the color 
factors, we found the algorithm to get a list of color factor (order r) from those of (order s) 
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and "T" (order r-s), {xo,A,B,C,D) < — (x[),A',B',C',D') x {x'i^,A" ,B" ,C" ,D"), 

xo < — x'o, (3.5) 

Ai — A'+A", (3.6) 

B< — B' + B" + AtA", (3.7) 

C ^ ( C[i), • ■ • ,C|,) , C[;) +2g^ + AfDj^), • • • + + A?D;;_,) ), (3.8) 

.s elements r-.v elements 



r elements 

D ^ (D'(i),...,D'(,),D'('i)+2A',...,D'(;_,)+2A'), (3.9) 

.V elements r-.v elements 

r elements 

where = Xq — Xq. It turns out that the resulting A and B remain single component integer. On the 
other hand, the resulting C and D are given by combinations of single prime and double prime ob- 
jects with some additional terms. Since A,B,C,D are all integer value and this operation is simple, 
the algorithm is suited for a symbolic operation and easily implemented in python script language. 
As a new ingredient in the implementation, we have to add the new components xo,A,5,C,Z) to the 
earlier list structure, 

lP = {al,---,a[,xo,Ak,Bk,Ck,Dk;fk). (3.10) 

Even in the SF with non-trivial color factor, the algorithm maintains a closed multiplication struc- 
ture, therefore it is applicable for any parallel transporter. 

To confirm and check the algorithm, we perform a one-loop calculation of the SF coupling. 
By calculating this quantity, we can check the two-point vertex, that is, inverse propagator. We 
investigate not only the plaquette gauge action considered before [||], but also the improved gauge 
actions including the rectangular loop [Q] and get consistent results. Furthermore, we compared 
with the hand derived three point vertex of the plaquette gauge action, available from a private note 
of Peter Weisz and confirmed consistency. To get further confidence in the implementation of our 
algorithm, we have to check the four-point vertex. To do so we need to do a two-loop calculation 
of the SF coupling [||, ^, and this will be reported in the future. 



4. Application I: L=T ±a lattice 

As a simple novel application of our algorithm, we perform a one-loop computation of the SF 
coupling on lattices with L = T + sa with 5 = ±1. Such lattices are motivated by considering the 
SF with staggered fermions. Due to the Dirichlet boundary in the time direction, we have to set 
T to be odd. Here, we will discuss only the gauge part on the lattice, not the staggered fermion 
part. More details about the latter are given in the contribution of P. Perez Rubio and S. Sint to this 
conference. 

When one takes the continuum limit in the standard SF, one sets the T /L = 1. The fact that we 
here have to set (r + sa)/L = 1 when taking the continuum limit in the tree-level 0{a) improved 
theory causes some change to the solution of the equation of motion, that is, the background field. It 
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Figure 1: One-loop relative deviation as a function of a/L. The left(right) panel shows tree(one-loop) level 
0{a) improved results 



is modified by lattice artifacts as compared to the standard one, and has to be extracted numerically 
from the equation of motion. The time dependence of the background field turns out not to be 
strictly linear anymore. Therefore we have to extend our algorithm to apply for an arbitrary time 
dependent phase 0(xo). For this purpose, we have derived another variant algorithm for the color 
factor which is similar to eq.(^^, and it will be shown in more detail in a future publication. 

In Figure we show the resulting relative deviation of the step scaling function to one-loop 
order, (a/L), given by 

5^u,a/L) = ^("'^/-^)-^(") =5i'\a/L)u + 0{u^), (4.1) 
a{u) 

a{u)=f{lL), u = fiL), (4.2) 

where k = 0(k = 1) is the tree (one-loop)level 0(a) improved case. In the plots, for comparison, 
we also include those of 5 = where L and T are the same. As a result, we observe that all three 
cases have similar absolute size of the lattice artifacts. 



5. Application II: A parameter for improved gauge actions 

As a next application, we apply the algorithm to the improved gauge actions including six 
link loops (not only for the rectangular type but also for the chair and 3-dimensional type actions). 
The loops are shown in Figure ^ with weight factors, cq, ci etc. and the weights are normalized 
by Co + 8ci + 16c2 + 8C3 = 1. We perform the one-loop computation of the SF coupling for the 
improved gauge actions and extract information about a ratio of lambda parameters for the pure 
SU(3) gauge theory between the lattice scheme (the various gauge actions) and the SF scheme, 
ALat/AsF- By combining the result of Asf/A^ in [|6|] for A/^ = 3 and N{ = 0, we summarize the 
values of ALat/Aj^ = ALat/AsF • Asf/A^^ in Table |l|. We observe rough consistency with old 
results. 

We thank the Deutsche Forschungsgemeinschaft (DFG) for support in the framework of SFB 
Transregio 9. 
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Co Cl C2 C3 

Figure 2: Loops with weight factors co,i,2,3 which contribute to the improved gauge action. 



Action 




C2 


C3 


ALat/A^s 


Plaquette 











0.0347109675... % 


Wilson RG 


-0.252 





-0.170 


2.34086(2) 


Iwasaki 


-0.331 








2.12455(5) 


DBW2 


-1.40686 








44.21(2) 


Symanzik 


-1/12 








0.1836938(4) 


Symanzik II 


-1/12 


1/16 


-1/16 


0.1782883(4) 


no name 


-1/12 


-0.1 


0.1 


0.1673674(8) 



Table 1: We show our results of ALat/A^^g in SU(3) gauge theory for various gauge actions, except for the 
plaquette gauge action ref. P] which is given here for completeness. Our values are roughly consistent with 
old results [|, |, [lO[ |ll], |l^ y], but in ref.[|l3| the authors do not show error estimates therefore, it is hard 
to compare. After the completion of this work, new results for several gauge actions have appeared in ref. 
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